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^ . Abstract 

> 

O 

(^ . Assuming three flavour neutrino mixing takes place in vacuum, we investi- 

^^ . gate the possibihty that the solar fe take part in MSW transitions in the 

^D ■ Sun due to A?7i|^ ~ (10""^ — 10~^) eV'^, followed by long wave length vac- 

~ uuni osciuauions on une way lo me izjann, Liiggereu uy la//12i {<-»: ^-^'"'32) 

^ . ~ (10~^^ — 10"^'') eV'^, ^'^ii S'lid Attt-Ii (Amgg) being the corresponding 

^^. neutrino mass squared differences. The solar u^ survival probability is shown 

^H. to be described in this case by a simple analytic expression. Depending on 

,^ , whether the vacuum oscillations are due to Ayrigi or Am'^2 there are two 

very different types of interplay between the MSW transitions and the vac- 
^ . uum oscillations of the solar f g . Performing an analysis of the most recently 

published solar neutrino data we have found several qualitatively new solu- 
tions of the solar neutrino problem of the hybrid MSW transitions + vacuum 
oscillations type. The solutions differ in the way the pp, Be and ^B neu- 
trino fluxes are affected by the transitions in the Sun and the oscillations in 
vacuum. The specific features of the new solutions are discussed. 
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1. Introduction 

It is well known that the solar neutrino problem [0-0] admits, in particular, two quite 
different neutrino physics solutions: one based on the old idea of Pontecorvo 0] that the 
solar Vf. take part in vacuum oscillations when they travel from the Sun to the Earth, and 



a second based on the more recent hypothesis @,|T0| of the solar i/g undergoing matter- 
enhanced (MSW) transitions into neutrinos of a different type when they propagate from 
the central part to the surface of the Sun. Both these solutions and the possibilities to test 
them and to distinguish between them in the future solar neutrino experiments, have been 



extensively studied in the simplest case of two- neutrino mixing (see, e.g., refs. |11-16| and the 
articles quoted therein). The vacuum oscillations and the MSW transitions of the solar v^. are 
characterized by the same two parameters in this case P|-pil|,P^: sin^ 2^, d being the neutrino 
(lepton) mixing angle in vacuum, and Am^ = m^ - vn\ > 0, where mi^2 are the masses of 
two neutrinos z/i 2 with definite mass in vacuum. The vacuum oscillation interpretation of 
the solar neutrino problem requires that ||16| 



5.0 X lO^^^eV^ ^ Am^ ^ 1.1 x 10"^°e1/^ (la) 

0.67 ^ sin'^29 < 1.0, (16) 



while the MSW solution is possible for |T6| 



3.6 X lO^^el/^ ^ Am^ ^ 9.8 x lO^^el/^ (2a) 
4.5 X 10"^ ^ sin'^29 ^ 1.3 x 10"^ (26) 

and for 

5.7 X 10~^e\/^ ^ Am^ ^ 9.5 x lO'^el/^ (3a) 

0.51 ^ sin'^2e ^ 0.92, (36) 

if the MSW transitions are into an active neutrino, ^^(t-) Q. The Am'^ solution interval (la) 



^For a detailed study of the MSW solution of the solar neutrino problem with solar Ug transitions 
into a sterile neutrino, Ug, see ref. ||lq| . 
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corresponds to values of the solar Ve oscillation length in vacuum comparable to the Sun- 
Earth distance. 

Detailed studies of the vacuum oscillation or MSW transition solution of the solar neu- 
trino problem under the more natural assumption of three flavour neutrino mixing are still 
lacking, partly because of the relatively large number of parameters involved. The three 
flavour neutrino mixing hypothesis implies the existence of three neutrinos Vk with definite 
vacuum masses, irik., k = 1,2, 3, and the relevant vacuum oscillation and/or MSW transition 
probabilities depend now on two different neutrino mass squared differences 0, Am2i > 
and Arrig]^ > 0, and on at least two mixing angles, say, 612 and ^13. The general features of 
the three-neutrino MSW transitions of solar neutrinos have been analyzed by many authors 
(see, e.g., refs. [p!9| , p0| ). Certain aspects of the possible three-neutrino MSW solution of the 
solar neutrino problem have been discussed as well |21|j2^ . However, comprehensive results 
based on the more recent and more precise data from the four solar neutrino experiments 
have been obtained only in specific cases when, for instance, Amji has a value in one of the 
solution intervals (2a), (3a), while AttIj^ >> 10^^ eV"^ E^HM. Under these conditions the 



solar z/e MSW transitions are essentially of the two- neutrino mixing type |p!U| , pU[ , although 
the regions of solutions can differ somewhat from those given in eqs. (2a) - (3b) [^,^. 

In the present article we study the qualitatively new possibility of a hybrid vacuum 
oscillation and MSW transition solution of the solar neutrino problem f\ This possibility is 
rather natural if three-neutrino mixing takes place in vacuum. 

One would guess a priori that such a solution would require Amli (or Am|2) and Am|^ 



^One can always choose mi < 1112 < m^ without loss of generality and we will work with this 
convention in what follows. 

^An earlier rather brief and qualitative discussion of combined solar i/g MSW transitions and long 
wave length vacuum oscillations in connection with the solar neutrino problem can be found in ref. 
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to have values in the regions specified by eq. (la) and by eqs. (2a) or (3a), respectively. 
This corresponds to the case 

Anil^ < Anil^ ^ Am^^, (4) 

or to 

Am^2 < ^"^31 - Am^p (5) 

Several patterns of neutrino masses can lead to the inequality (4) or (5) and to the 
requisite values of Ar/igi (o^' ^^2) ^'^d ^"^ii- Eqs. (la), (2a) or (3a) and (4) (or (5)) are 
compatible, for instance, i) with the existence of three quasi-degenerate neutrinos, mi = 
7712 — "^3; with a very different splitting between their masses, {1712 — mi) <^ (777,3 ~ "^i); ii) 
with a neutrino mass spectrum consisting of two close masses which are very different from 
the third one, mi = 777,2 ^ ^^3, and iii) with a hierarchical neutrino mass spectrum, r7ii -C 
7772 ^ ^^3- In the first case neutrinos can have masses in the cosmologically relevant region, 
mi ~ (1 - 2) eV, with (mg - mi) ~ (10"^ - 10"^) eV and (ma - mi) ~ (IQ-^^ _ xg-io) eV, 
while in the other two cases the neutrino masses are too small to be of significance for 
the solution of the dark matter problem: ma ~ 10^^ eV and 7772 ~ 7r7i ~ (7772 — ttii) ~ 
(10^^ — 10"^) eV (case ii) or 7712 ~ (10^^ — 10^^) eV (case iii). In case ii) the pattern 
777i ^ 7r72 = 7r73 with 7712,3 ~ 10~^ cV and (7723 — 7772) ~ (10^^ — 10^^) eV is also possible. It 
corresponds to eq. (5) and arises in the Zee model of neutrino mass generation p6|-p8| and 
its various GUT generalizations. Obviously, if A77I21 (A7r7|2) and A777|^ lie in the intervals 
(la) and (2a) or (3a) respectively, there will be no observable neutrino oscillation effects 
due to Am,2i (A7r732) and A7773^ in the oscillation experiments performed with terrestrial 
(reactor, accelerator, etc.) or atmospheric neutrinos. 



"^If the massive neutrinos fj are Major ana particles, the neutrinoless double beta decay would be 
allowed and, for 77ij ~ (1 — 2) eV, can proceed with a rate which is in the range of sensitivity of 
the new generation of experiments searching for this decay (see, e.g., [p5| ). 



All the different patterns of neutrino masses mentioned above can arise in gauge theories 
of electroweak interactions with massive neutrinos, and in particular, in GUT theories and 
in the superstring inspired models (see, e.g., p9|). 

This paper is organized as follows. Section 2 is devoted to the discussion of the average 
solar z/e survival probability, -P(z/e -^ z/g), when three-neutrino mixing takes place in vacuum 
and the solar neutrinos undergo MSW transitions in the Sun, which are followed by long 
wave length neutrino oscillations in vacuum. It is shown, in particular, that the probability 
Pl^e — > i^e) is described by a simple analytic expression and that inequalities (4) and (5) 
imply two very different types of interplay between the MSW transitions and the vacuum 
oscillations of the solar z/g . In Section 3 the hypothesis that the solar z/g take part in MSW 
transitions and long wave length vacuum oscillations is confronted with the solar neutrino 
data. Several new hybrid MSW transition + vacuum oscillation (MSW + VO) solutions 
of the solar neutrino problem are found and their specific features are discussed. Section 4 
contains our conclusions. A derivation of the analytic expression for the probability -P(z^e -^ 
Ve) used in the present study is given in the Appendix. 

2. The Solar z/g Survival Probability 

It is possible to obtain a very simple analytic expression for the solar z/g survival probabil- 
ity in the case of interest, namely, when the solar z/g undergo matter-enhanced (MSW) transi- 
tions in the Sun and vacuum oscillations on their way from the surface of the Sun to the Earth. 
The vacuum oscillations are supposed to proceed for neutrinos having energy E' ~ 1 MeV 
with an oscillation length comparable to the Sun-Earth distance, Rq = 1.4966 x 10® km. 

We assume three-flavour neutrino mixing takes place in vacuum, 
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Wi> =J2U*kWk>, / = e,/i,r, (6) 

fe=i 

where {ui > is the state vector of the (left-handed) flavour neutrino ui having momentum 

jf, |z/fc > is the state vector of a neutrino u^. possessing a definite mass rrik and momentum 

jf, irik 7^ TTLj, k ^ j = 1,2,3, nil < "^2 < ""^3, and U is a 3 x 3 unitary matrix - the 



lepton mixing matrix. The neutrinos Uk are assumed to be stable and relativistic. We shall 
suppose also that the matter-enhanced transitions of the solar z/g into z/^ and/or u^ in the 
Sun are associated with Am'^i, where, as usual, Am?- = mf — m'j, while Am'^i (or A'm'^2) 
is responsible for the long wave length (~ 1.5 x 10® km) Ue ^->- ly^^r oscillations taking place 
between the Sun and the Earth. This implies that either relation (4) or (5) is valid . We 
shall consider in the present analysis a somewhat wider range of values of Ar7i2i {^''^32) ^^^ 
Amg^ than is given by the two- neutrino mixing solution intervals eqs. (la) and (2a) or (3a), 
respectively: 

10"^^ eV"^ < Aml^ ^ 5.0 x 10"^° eV\ (7) 

10"^ eV^ < Am^i < 10"^ eV\ (8) 

Under the above conditions and assuming the inequality (4), Am'^i ^ Am|]^, holds, the 
analytic expression for the solar z/g survival probability of interest, -P(t'e -^ ^e)-, can be easily 
deduced from the general expression for P{U(, -^ u^) derived in ref. ||2^ (see eq. (29) in pO| ) 



for the relevant case of three (flavour) neutrino mixing. It has the following simple form: 



/Q1 \ 

P(z/e -^ Pe, tE, to) = P2MSW\^e -^ ^e\ ^0, ^o) 

+ ]~f\jf\2 \^-P?yo{^^ ^ ye\tEM\ W,^f-P^2^^sw{^e - Ve.t^M)\- (9) 
1 - l\Ue'iY 

Here -P(z/e -^ Ve'-, ^e, to) is the average probability that a solar u^ having energy E' ~ |"p^| = p 
will not be converted into t'^ and/or Ur when it propagates from the central part of the 
Sun, where it was produced at time to, to the surface of the Earth reached at time tg, 
P2Mswi^e — * VeitQ,to) is the average solar z/g survival probability when the z/g undergoes 
two-neutrino MSW transitions due to Arrigj^ and f/es 7^ on its way to the surface of the 
Sun, t0 being the time at which the z/g reaches the surface of the Sun, and P2voi^e —>■ 
^e'ltEito) is the solar z/g survival probability for a z/g taking part in two-neutrino vacuum 
oscillations due to Am\^ and f/gi 7^ 0, Ue2 7^ 0, on its way from the surface of the Sun to the 
surface of the Earth. The two-neutrino MSW transition and vacuum oscillation probabilities 
-^2M5'H/('^e — >■ i^e',tQ,to) and -R^vm^e — ^ ^'e; ^E, ^0) are given by the well-known expressions 



liPiO^: 



,1 
2 ■ '2 



PiMSwi^e^l^e^t^M ^ P^'^sW = ' + {^ - P'^u)) ^Os2e^, COs2eTs{to) , (10) 



Pivoi'^e ^ 'ye;tE,te) = P^S = 1 " ^ siii^ 2^i2 (1 - cos27rA). (11) 

In eqs. (10) and (11) PL-^^^ is the so-called "jump" (or "level-crossing") probability, 



COs2^13 = l-2|f/e3|'>0, sin2^i3 = 2|f/e3|Vl-|f^e3P > 0, (12) 

6^ is the two-neutrino mixing angle in matter, which in vacuum coincides with 6*13, 

^(l-^)2 + tan2 2^i3 

where A^e(^o) is the electron number density at the point of u^ production in the Sun and 
j\jres _ /S.ml'^ cos29i3/{2p\/2Gf) is the resonance density, 

R = [Ie — to) is the distance traveled by the neutrinos in vacuum, R = Rq, and L21 = 
47rp/A?Ti2i is the oscillation length in vacuum associated with Arri^i- 
Several comments are in order. The requirement 

\Ue3\^ < 0.5 (15) 

in eq. (12) ensures the possibility of a resonant enhancement of the solar z/g transitions in 
the Sun |4|,|. 

For the values of A?7i|^ of interest, eq. (8), the oscillating terms in the probability 



P(z/e -^ Ue] tE, to), associatcd with the solar z/g MSW transitions in the Sun ||35| are rendered 
negligible [^ by the various averagings one has to perform when calculating the effects of 
the transitions and oscillations on the signals in the solar neutrino detectors. 

Further, the probability -P(t'e -^ ^e',tE,to), eq. (9), depends only on the absolute values 
of the elements of the lepton mixing matrix U, forming the first row of U, more precisely, 
on |f/efcP; k = 1,2,3. Of these only two are independent since the unitarity of U implies: 
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|[/eip + |f/e2p + \Uez\'^ = 1- It follows froiii eqs. (12), (14) and the preceding remark that 6i^ 
and 6i2 are independent parameters. 

We note next that at least three different expressions for the "jump" probabihty -P/31) 
entering into formula (10) for the MSW probability have been proposed [^,0. We will use 
the one derived in ref. [^ in the exponential density approximation: 

g-27rro^^sin^6»i3 _ g-27rro^^ 

where tq is the scale-height of the variation of the electron number density Ne along the 
neutrino path in the Sun. Unlike the expressions proposed in ref. [^, the one given by eq. 



(15) describes correctly the (strongly) nonadiabatic MSW transitions of solar neutrinos for 
values of sin^ 2^13 ^ (0.2 - 0.3) |31]j33,|33. 

Let us comment on how eq. (9) with probabilities P2MSW ^^"^ P-lvo given by eqs. (10) 
and (11) can be obtained from eq. (29) in ref. [^. The analysis leading to eq. (29) in |^ 
is valid, in particular, under the condition Am'^i ^ Aml^. This condition is fulfilled in the 
case presently considered, which is specified by eqs. (4), (7) and (8). We observe first that in 
this case one has, in the notations of ref. [^], £12 = 0, ^Pisit'o) — (and therefore £13 = C13, 
Si3 = S13) and S23 = in eq. (29); as a consequence, the latter simplifies considerably. We 
notice next that for Am2i having values in the interval (7) the times t and tg in eq. (29) in 
pO| can be considered to be respectively the times Ie and t© with the meaning they have in 
eq. (9) above: the derivation of eq. (29) in |^ for -P(z/e —>■ Uf.) is valid in this case as well. 
Finally, i) the angles ipi2 and ipis in eq. (29) in |20[ coincide with the angles 612 and ^13 
defined by eqs. (12) and (14) above, ii) the probabilities P^Ug — > Ug] t, to), Pni^e —^ ^e', t'o, to) 
and Piii^e -^ ^e]t,tQ) in eq. (29) coincide for t = Ie and tg = t© with the probabilities 
P{ue -> i^e;tE,to), P-^Mswi^e ^ z/e;t0,to) and P-^voi^e -^ Ve^tE^to), respectively, and iii) 



the term Re [i?//(tQ, to)] in eq. (29) in the case of interest is given, as it is not difficult to 
show using the results in refs. |p2| , |35| , ^ , by the following expression: 



Re [Rnit'^M] = -{\- ^31)) sin 2^13 cos 2^r3(^o). (16) 



Given the above observations it is easy to obtain expression (9) for P(t'e -^ i^e',tE,to) from 
eq. (29) in I^O] [|. A proof of eq. (9) which does not rely on eq. (29) in |^0[ and is based 



only on some of the intermediate results of the analysis performed in |2^ is given in the 
Appendix. 

Using eqs. (10), (11) and (12) we can rewrite the second term in eq. (9) in a form which 
proves convenient for our later discussion: 

P{u, ^ u,;tE,to) = P^mIw-I sin22^i2 (l-cos27r-^) cos2^i3[^ + (i-i^3i)) cos 20^3(^0)] 

(17) 

Before proceeding further let us point out that since P2MSW — l^esP = sin^ 6'i3 §] and 
|f^e3p < 0.5, for fixed values of the MSW transition parameters one always has 

P(z/e - ^e; tE, to) > ^^^^^^^2 i^ ' l^^sl' ' HmIw] = Sm' ^13[^ " (^ " ^31)) COS 20^(^0)] 

(18) 

as it follows from eq. (9). This lower limit constrains effectively -P(t'e ^ i^e]tE,to) only for 
sufficiently large values of sin^6'i3. Setting P2voi^e. -^ ^^e'ltE-ito) = 1 in eq. (9) we get an 
absolute upper limit on P{i'e — *■ i^e'-,tE,to) for given P2MSW ^^cl |f/e3p: 

P(z/e -> Z/e; tE, to) < HMlwi'^e ^ J^e', t©, to). (19) 

The three-neutrino mixing solar z/g survival probability in the case of interest is always 
smaller than the corresponding two-neutrino MSW probability. 

Eqs. (18) and (19) determine the range in which the probability -P(z/e -^ z/e;t£;,to) can 
oscillate as a function of R/p due to the presence of the vacuum oscillation component P2V0 
in it. Let us consider several specific cases relevant to our further discussion. 



^Note that there is a misprint in eq. (29) in [^]: the probability Piii^e -^ ^'e'^i^t'o) in the 



expression in the curly brackets multiplying the term Re [i?//(iQ,io)] in the third row should read 

PliVe -^ l^fi] t, t'o) = 1- Plil'e -^ ^e] t, tg)- 
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If the solar v^ undergo extreme nonadiabatic MSW transitions in the Sun (p/Ar/ig]^ is 
rge"), one has [0 P^^^ = cos^^ia, P?Ssw = 1 - | sin' 2^13 and 

U,,\^ = Sin^^i3 < Pi'^e ^ iye;tE,to) < \U,s\' + (1 " It/esH' = sin^ ^13 + COS^^ig. (20) 



We get the same result, eq. (20), if A^e(^o)(l — tan26'i3) ^ < N^^^, i.e., when p/Aml^ is 
sufficiently small. In this case |3^ the z/g MSW transitions are adiabatic (-P^i) — 0), the 



Ue oscillate in the Sun due to Aml^ as in vacuum: cos26'5^(to) — cos 26*13 and P2MSW — 
1 — |sin^20i3. If, however, in the "small" p/Arn^^ region the vacuum oscillation term 
cos27ri?/L2i averages out so that effectively -P2V0 ~ -'- ~ l/^ sin^ 26^12, one obtains: 

P{U, -^ Ue]tE,to) = \Uel\^ + \Ue2\' + Wes]' = 1 - ^ siu^ 2^i2 COS^ ^13 - ^Sin2 2^i3, (21) 

which is the average three-neutrino vacuum oscillation probability. 

If the solar Ue MSW transitions are adiabatic (i.e., -P/31) — 0) and cos20']^(to) — —1 
(i.e., Ne{to)/N^'^^ ^ 1, tan 2^13, solar neutrinos are born "above" and "far" (in Ne) from the 
resonance region), one has P2MSW — l^esP and therefore 

PiUe -^ l^e\ tE, to) = P^mIw = Pes]^ = Sm' ^13- (22) 

Under the above conditions the z/g state in matter at the point of u^. production (in the Sun) 
essentially coincides with the heaviest of the three neutrino matter-eigenstates, which con- 
tinuously evolves (as the neutrino propagates towards the surface of the Sun) into the mass 
(energy) eigenstate |z/3 > at the surface of the Sun. As a consequence, vacuum oscillations 
do not take place between the Sun and the Earth and P{h'e -^ ^e',tE,to) coincides with the 
probability to find u^. in the state lz/3 >. 

Let us discuss next the properties of the probability -P(z/e -^ ^e]tE,tQ) if inequality (5), 
i.e., a Zee model type relation between Am'^2 and Am"^^ (Ar/igi)) is valid and Am'^2 is assumed 
to give rise to solar neutrino long wave length vacuum oscillations on the path between the 
Sun and the Earth. The probability -P(z/e -^ Ue] tE, to) in this case has the form: 

P (t'e — >■ J^e] tE, to) = P2MSwi^e ~^ ^e'j ^0? ^o) 



+ ^JZf\2 [l-^2yo(^e ^ ^e.tE.tr^)] [W el? - P^ISW i^ e ^ ^e.tc.t,)] (23a) 
i Z|Cei I 

ASh^-^ sin22^23 (l-cos2vrA) siri^ e[,[]^- {]^- P[,^))cos2d[^{to)]. (236) 



The angles ^23 and 6*^3 in eq. (23b) are determined by 

We?.? \Ue2\ 
\\Ue,? + \Ue2?? 

and 



\TJ |2 irr |2 
sm 26'23 = 4 ,|.., |„ , ,■■, |^.„ , (24) 



cos2e;3 = 2|f/eir-l >0, sin2^;3 = -2|t/ei|Vl- If^eiP <0, (25) 

and L32 = 47rp/Am32; the requirement 

pel? > 0.5 (26) 

guarantees the possibihty of resonant enhancement of the MSW transitions in the Sun. 
The two-neutrino mixing MSW probabihty P2Mswi^e ~^ ^e',tQ,tQ) = P2MSW i^ given by 
eq. (10) in which the angles 6^13, 9^{tQ) and the "jump" probability Pki) are replaced 
respectively by 9[^, O^ito) and -P/13), where O^^ito) and P'n'j,-) are defined by eq. (13) and by 
eq. (15) in which ^13 is substituted by 9'^^. The two-neutrino vacuum oscillation probability 
Avm^e — > ^^e'ltEito) = P2V0 ^^^^ ^e obtained by replacing in eq. (11) 612 and L21 with 6*23 
and L32. 

Expressions (23) follow Q from eqs. (9) and (17) if we interchange the indices 1 and 3 
in the quantities \Uej? and Am^^, or equivalently, if we replace |f/e3(i)P with |f/ei(3)P and 
Am-g-^ (Ar?22i) with Am^g (Ar?7,23) in all the terms entering into eqs. (9) and (17), except in 
the "jump" probability -P/31), and take into account the definitions of 61^ (eq. (12)) and ^^3 
(eq. (24)). It can be shown that P'{i2,) coincides with the expression for -P/31), eq. (15), in 
which ^13 is changed to O'-^^,- 



^In the particular case of the Zee model lepton mixing matrix U |2^,27| in which two of the three 
mixing angles have specific values and all elements of U are real, expression (23a) for P^ {I'e -^ 



V(,]tE,tQ) reduces to the one derived in ref. |28| in the indicated case. 
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Although expressions (9) ((17)) and (23a) ((23b)) formally look quite similar, the two 
probabilities P{i'e ~^ J^ei^c^o) aiid P^{ve —* i^e]tE,to) actually have very different proper- 
ties. 

Expression (17) implies that if Am'^i <^ Am|]^, a genuinely hybrid small mixing angle, 
sin^ 6'i3 ^ 1, MSW transition + vacuum oscillation solution(s) of the solar neutrino problem 
is (are) in principle possible. In contrast, if sin^ 9[.^ <^ 1 the contribution of the vacuum 
oscillation probability -P2V0 i^ P^i^e —* i^e'-,tE,to) is strongly suppressed and P^ii^e — ^ 
Ue'jtEyto) = P2Mswi^e — ^ ^e',tQ,to). Thus, if inequality (5) holds, hybrid MSW + VO 
solutions can exist only for large values of the MSW mixing angle ^^3. 

Similarly to -P(z^e -^ ^e',tE,to), the probability P^iye — ^ ^e',tE,to) is limited from above 
by P2MSW- P^i^e — ^ i^e',tE,to) < P2Mswi^e — ^ i^e',tQ,to) f\. This follows from the fact that 



0.5 < \Uei\^ < 1.0 and [|t/eip- Amsh/] > 0- However, the lower limit on P^(i/e -^ i^e] ts, t 
is physically very different from the lower limit on P{ve — ^ '^e',tE,to), eq. (18): 



0) 



(27) 



P'{U, -. U^-tE.to) > J, ?J , [PSlw-i^-lUell')] = COs2^;3[- + (--P('i3))cOs2^;™(to)] 



In the two "asymptotic" regions of "small" (adiabatic MSW transitions) and "large" (ex- 
treme nonadiabatic MSW transitions) values of p/Ar/il^ the probability P^{ve -^ i^e'itE^to) 
has the form: 

P^(z/e ^ ^e;tEM) = \Uel\^ + (1 " Wel^fP^^k^e "^ ^^e;^^,^©). (28) 

From eq. (28) we get the following band of allowed values of -P^(t'e -^ ^e]tE,to), which is 



For the specific case of the Zee model mixing matrix U |26,27| this upper bound was obtained 



inref. [g8|. 



^We have used the upper bound P2MSPV ^ ^ " 2\Uei\'^{l - |f/eiP) = 1-1/2 sin^ 261^3 to get the 
last inequality. 
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the same for the two regions: 

\Uel\^ = COS^^^g < P^(z/e ^ iy,;tE,to) < \Uelf + (1 - \Uel\'')^ = COS^ ^^3 + sill^ ^^3 . (29) 

The corresponding band for the probabihty -P(z/e — > z/g; tE, to) is given in (20). 

Finally, if the MSW transitions in the Sun are adiabatic, PA3-) = 0, and cos 29^(to) = —1 
(HmIw = sin^^i3), one has: 

< P^iUe -^ Ve] tE, U) = (1 " \U ex\^)P?vl[Ve -^ Ve] tE, t©) < sin^ ^^3 , (30) 

to be compared with eq. (22). 

It follows from the above discussion that if Ar?7,2i ^ ^"^su the amplitude of the oscil- 
lations of -P(z/e -^ ^e'ltE-ito) (due to the -P2V0 term) as R/p varies is equal in the regions 
of "large" and "small" p/ Arrag^ to cos^ ^13 sin^ 26*12 and can be maximal; for sin^ 2^12 ~ 1 
it is always greater than 1/4. At the same time, this amplitude is strongly suppressed and 
vacuum oscillations practically do not take place (even if sin^ 2^12 ~ 1) in the region of 
the adiabatic MSW transitions where cos 2^]^ (to) = —1 and P2MSW ^^^ its minimal value, 
P2MSW ~ sin^^i3. In contrast, if the Zee model type relation Ar?7,|2 <^ Arrig^ holds, the 
vacuum oscillations are suppressed in the two "asymptotic" regions: the amplitude of the 
oscillations is equal to sin^ 0'^^ siv? 2^23 and is always constrained to be smaller than 1/4. In 
the region of adiabatic MSW transitions, where cos 26*1^(^0) — —1 and P2MSW ~ ^™^ ^i3' the 
probability P^ii^e -^ i^e',tE,to) oscillates (due to -P2V0) ^^ a function of R/p with an ampli- 
tude sin^ 9[^ sin^ 26*23, which does not exceed 0.5 sin^ 26^23- Obviously, the vacuum oscillations 
can be important in this case only for sufficiently large values of sin^ ^'^3. 

Because of the aforementioned differences in the properties of the probabilities -P(t'e ~^ 
^e] tE, to) and P^{ve -^ ^e, tE, to), the case (4) of relations between the neutrino mass squared 
differences provides a much richer spectrum of possibilities of genuinely new solutions of the 



^Note that the band of allowed values of P'^{ve -^ ^e',tE,to) in the "large" p/Atti^i region shown 



in Fig. 2 in ref. |28| does not correspond to eqs. (28) and (29). 
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solar neutrino problem of the hybrid MSW transition + vacuum oscillation type, than the 
case of relations (5). 

3. Hybrid MSW Transition + Vacuum Oscillation Solutions 
of the Solar Neutrino Problem 

We have analyzed the most recently published data from the four solar neutrino ex- 
periments Il|-3 searching for solutions of the solar neutrino problem of the hybrid MSW 
transitions + vacuum oscillations type. Only the case of inequality (4) with Amg^^ and 
Am2i having values in the intervals (7) and (8) was studied. 

The analysis we performed is based on the analytic expression eq. (9) for the average 
solar i^e survival probability, -P(z/e — > Ve]tE-ito)-i with P2Mswi^e — > '^e;^0,^o) given by eqs. 
(10), (13) and (15). As is well known, the solar neutrino flux consists of several components, 
six of which are relevant for the interpretation of the results of the solar neutrino experiments 
(see also, e.g., 0,0]): the pp, pep, ^Be, ^B and the CNO (two components). We utilized 
the predictions of the solar model of Bahcall and Pinsonneault from 1995 with heavy element 
diffusion |3^ for the pp, pep, etc. neutrino fluxes in this study. The estimated uncertainties in 



the theoretical predictions for the indicated fluxes |3^ were not taken into account. For each 
given neutrino flux the probabihty P{ye -^ ^e',tE,to) was averaged over the corresponding 
region of u^ production in the Sun. Since the data analyzed was accumulated over a period 
of several years, P{i'e -^ ^e',tE,to) was also averaged over an interval of time equal to one 
year [|ll|]. The probability -P(z/e -^ ^e',tE,to) depends on the time of the year ty through 
the dependence of P2voi^e "^ ^e]tE,tQ) on the distance between the Sun and the Earth, 
ts — tQ = R, which is a function of ty-. 

R = R{ty) = i?o [1 - e cos 27r^], (31) 

where e = 0.0167 is the ellipticity of the Earth orbit around the Sun and T = 365 days. 
The solar neutrino fluxes are functions of R~^{ty) and their dependence on ty was effectively 
accounted for when we performed the time averaging of the probability P{ve -^ '^e', tE, to). 
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The results we shall report have been obtained utilizing the x^~ method. The regions 
of the relevant MSW transition and vacuum oscillation parameters allowed by the data 
were determined by requiring that x^ < 3.841. This choice is motivated by the following 
arguments. In the x^~analysis in the case under study there are zero degrees of freedom: we 
have four experimental results and the theory tested contains four parameters (6'i2, Atti^i, ^13 
and Ar/ig]^). Nevertheless, the existence of regions of values of the parameters for which the 
X^— function has a distinct minimum and a sufficiently low value at the minimum is a strong 
and unambiguous indication that the mechanism of solar z/g flux depletion considered gives 
a good quality of the flt of the data and therefore provides a solution of the solar neutrino 
problem. At the same time, in the case of zero degrees of freedom it is impossible to assign 
a precise confldence level to the x^ = 3.841 contours. It should be noted, however, that in 
grand unifled theories with massive neutrinos the quantities Am2i and Am|^ are typically 
not independent: their ratio is a function of the ratio of quark or charged lepton masses 
Q. If such a relation holds there would be three independent parameters in the theory, and 
correspondingly one degree of freedom in the x^— analysis. Under these conditions x^ = 3.841 
corresponds to 95% C.L. 

We have searched for MSW + VO solutions of the solar neutrino problem for values of 
sin^ 2612 and sin^ 2^13 from the intervals: 

0.20 ^ sin2 2^i2 < 1.0, (32) 

10~^<sin2 2ei3 < 1.0. (33) 

The intervals (32) and (33) are considerably wider than the corresponding ones in the case 
of purely two-neutrino mixing (2z/) vacuum oscillation and MSW transition solutions (see, 
e.g., [p!^ -[l6[|). We did not perform a complete scan of the 4-parameter space deflned by eqs. 
(7), (8), (32) and (33). We were primarily interested in solutions which do not converge 



^'^In certain classes of theories some of the lepton mixing angles are also functions of the charged 
lepton or quark masses. 
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continuously (via regions where x^ < 3.841) to the 2v vacuum oscillation or MSW transition 
solution in the limit of sin^ 2^13 -^ or, correspondingly, of sin^ 2^12 ^ or Am^g ~^ 0- 
Nevertheless, we have found five different regions of values of the four parameters, for which 
X^ < 3.841. These regions (solutions) are denoted by A, 5, C, D and E below. Sections 
of the indicated regions in the planes of the MSW transition and of the vacuum oscillation 
parameters Arrig-^ — sin^ 2^13 and Am'^i — sin^ 2^12 are shown respectively in Figs, la, lb and 
2a, 2b. Shown in Figs, la - 2b are also the points in which the x^~function has a minimal 
value. 

We shall discuss next the physical features of the solutions found. A detailed analysis of 
the implications of the MSW transition + vacuum oscillation solutions for the solar neutrino 
experiments Super-Kamiokande, SNO, ICARUS, BOREXINO, HELLAZ and HERON will 
be given elsewhere. Here we will limit ourself only to qualitative remarks. 

Solution A. For this solution (see Figs, la and 2a) the vacuum oscillation parameters 
lie in the intervals 

4 X 10-^2 eV^ ^ Amj-^ ^ 8 x 10^^^ eV^, (34a) 

0.65 ^ sin2 2^i2 < 1.0, (346) 

while the region formed by the values of the MSW transition parameters extends from the 
domain of the 2z/ nonadiabatic solution at Aml^ = 5.0 x 10^^ eV"^, sin^ 26^13 = 8 x 10~^, 
first to larger (smaller) values of Aml^ (sin^ 2^13) reaching Am'^i = 10^^ eV^ and sin^ 2^13 = 

3.1 X 10"^, and then to larger values of sin^2^i3 up to 0.5, spanning approximately 3 orders 
of magnitude in sin^ 2^13 at practically constant Aml^ = (1.1 — 1.3) x 10^^ eV'^. The 
minimum x^ value for this solution is Xmin — 0.071 and corresponds to (Am2i, sin^2^i2; 
Am^i, sin2 2^i3) = (5.6 x 10^^^ eV^, 0.98, 4.2 x 10"^ eV^, 10"^) (the black triangle-down in 
Fig. la). In the "horizontal" region where sin^ 26^13 ^ 10"^ and Am^^ = (1.1-1.3) x 10"'' eV^ 
there exists a local minimum with a somewhat larger value of Xmin'- one finds Xmin — l-O- 
This minimum is reached for (Am2i, sin^2^i2, Am|^, sin^2^i3) = (5.4 x 10~^^ eV'^, 1.0, 

1.2 X 10"^ eV^, 0.14) (the black dot in Fig. la). 
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For p = E > 5 MeV and Am^^ < 8.0 x IQ-^^ eV^ we have L^^ > 1.6 x 10^ km, 
cos{2-nR/Ll^) > 0.82 and, as follows from eq. (11), Pivoi'^e ^ Ve^tE.to) > 0.91. Thus, for 
®B neutrinos with energy E > 5 MeV one has in the case of solution A: -P(z/e —>■ Ve'-, ts, to) — 
P2Mswi^e -^ ^e] ^Q, to); i-G-, uiost of the ^B ueutriuos undergo only MSW transitions |^. The 
MSW transitions of the ^B neutrinos having energy E > 5 MeV are adiabatic for values of 
Aml-^ ^ (1.1 - 1.3) X 10-^ e1/2 and sin^ 2^13 ;^ (3.0 - 4.0) x 10"^ 13 from the "horizontal" 



region of the solution (see Fig. la). They are nonadiabatic for values of Am'^-j^ and sin^ 2^13 
from the remaining part of the allowed region. 

For the pp and the major part of the ^Be neutrinos we have E < 0.41 MeV and E = 
0.862 MeV, respectively, and for Ar77,|]^ ^ 1.3 x 10"^ eV"^, their energies fall in the "small" 
p/Am3^ domain (see Section 2) where 



P{iye -> Ue] tE, to) = sin'' 6113 + COS^ 613 Pi^oi'^e "^ J^e] tE, t©) 



(35) 



with sin^ ^13 ^ 0.15 for the solution under discussion. Thus, in contrast to the main fraction 
of ^B neutrinos, the pp and ^Be u^, do not undergo resonant MSW transitions but take 
part in vacuum oscillations between the Sun and the Earth. Actually, the ^Be neutrino 

(21) 

energy of 0.862 MeV is in the region of the first minimum of -P2V0 ^^ -^ decreases from the 

(21) 
"asymptotic" values at which P2V0 — -'-' while the interval of energies of the pp neutrinos, 

relevant for the current Ga-Ge and the presently discussed future solar neutrino experiments 

(HELLAZ, HERON), 0.22 MeV ^ E < 0.41 MeV, is in the region of the first maximum of 

P2V0 ^s E decreases further. The dependence of the probability P(z/e — > VeitE-,to) on the 

neutrino energy E in the case of solution A is illustrated in Figs. 3a - 3c. 

For 5.0 X 10~^ eV"^ < Amg^ ^ 1.3 x 10~^ eV"^, i.e., in the minor sub-region of solution 

A, which overlaps with the region of the 2v nonadiabatic solution (see Fig. la), the 0.862 

MeV ^Be neutrinos take part either in MSW transitions and vacuum oscillations or in MSW 



^^Obviously, if, for instance, Am^i = 4.0 x 10 ^^ eV'^ , the same result will be valid for the 
neutrinos with E > 2.5 MeV. 
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transitions only. This small sub-region converges continuously to the region of the 2z/ nona- 
diabatic solution when sin^ 9i2 -^ and we are not going to discuss it further. Note that the 
domain of solution A in the Am2i — sin^ 9i2 plane shown in Fig. 2a corresponds to values of 
Amli and sin^ ^13 which lie outside the sub-region in question. 

For solution A, the Kamiokande signal and the signal due to the pp (^B) neutrinos in 
the Ga-Ge (Cl-Ar) experiments are smaller approximately by factors of ~ (0.33 — 0.42) and 
~ (0.70 — 0.80) (~ 0.30), respectively, than the corresponding signals predicted in ref. [ ^ . 
The '''Be u^ flux is suppressed rather strongly - by a factor of ~ (0.13 — 0.14), with respect 



to the flux in the solar model of ref. ||39|| . 



Let us note that the solution A region in the ^tn^^ — siv? 612 plane of the vacuum oscil- 
lation parameters is quite similar to the region of the "low" '^B v^, flux 2v vacuum oscillation 
solution found in ref. [|ll|]. The latter is possible for values of the ^B neutrino flux which are 
lower by a factor of 0.35 to 0.43 (of 0.30 to 0.37) than the flux predicted in [^ (in [^ and 
used in the present study). Note, however, that the Xmm ^^^ ^^e indicated purely vacuum 
oscillation solution, xlnin — ^-4 (2 d.f.) |TT|, is considerably larger than the value of xlnin ^^^ 
solution A. The two solutions differ drastically in the way the ^B neutrino flux is affected by 
the transitions and/or the oscillations. 

The implications of solution A for the future solar neutrino experiments aimed at detec- 
tion and studies of the ^Be and/or pp neutrinos like BOREXINO, HELLAZ or HERON, are 



practically the same as those discussed in detail in ref. |jTT| for the case of the corresponding 
purely vacuum oscillation "low" ^B neutrino flux solution. However, the implications for the 
Super-Kamiokande, SNO and ICARUS experiments which are sensitive only to ^B neutrinos 
are completely different. They will be discussed elsewhere. Let us mention only here that 
in the case of solution A: i) the spectrum of ®B neutrinos will be strongly deformed (see 
Figs. 3a - 3c), ii) the magnitude of the day-night asymmetry in the signals in the indicated 
detectors can be very different from that predicted in the case of the 2z/ MSW solution (see, 
e.g., [^), and iii) the seasonal variation of the '^B u^ flux [|r^,|r^ practically coincides with 
the standard (geometrical) one of 6.68%. 
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Solution B. For solution B (see Figs, la, 2a and 2b) the allowed region in the Arn^i — 
sin^ 26'i3 plane of the MSW transition parameters, 

5.1 X 10"^ eV^ ^ Am^i ^ 1.2 x 10"^ eV^, (36a) 

3.2 X 10"^ ^ sin2 20i3 ^ 6.6 x 10~^ (366) 

looks approximately like the region of the 2z/ MSW nonadiabatic solution shifted as a whole 
to somewhat larger values of Arrig^ (by a factor ~ 1.3) and to smaller values of sin^ 2^13 (on 
average by a factor ~ 1.6). The two regions have a small common sub-region (see Fig. la) 
and our further discussion does not extend to this sub-region of overlapping. Solution B can 
be regarded as an "improved" MSW transitions -|- vacuum oscillations version of the purely 
2z/ MSW nonadiabatic solution. 

In the Amli— sin^ 2^12 plane the region of solution B extends in Am|^ from 1.8x10"^^ eV"^ 
at least up to 10~^ eV"^ (we did not perform a search for allowed values of Am^i beyond 
10^^ eV"^) and from 0.15 to practically 1.0 in sin^2^i2: 

0.15 ^ sin2 2^i2 ^0.98. (37) 

For Amli from the domain ~ (0.5 — 1.0) x 10^^° eV"^ of the 2z/ vacuum oscillation solution 



Tl|,|T6|, however, solution B takes place for values of sin 2^12 which are systematically smaller 



than the values of the same parameter in the 2z/ vacuum oscillation solution. As is seen in Fig. 
2b, in the Am^i — sin^ 26*12 plane there are two marginally disconnected regions, which we will 
call "down" and "up": the separation "line" is approximately at Am2i — 1-5 x 10"^" eV"^. 
Except in the region where Am^i = (1.0 — 2.0) x 10^^" eV"^, the vacuum oscillation mixing 
parameter is relatively small for solution B: one has sin^ 26^12 ^ 0.8, and for most of the 
allowed values of Am^i actually sin^ 26*12 ^0.7. 

The minimum x^ value for the solution under discussion is even smaller than for solution 
A: we have xLu = 0-025 at (Am^i, sin2 20i2, Am|i, sm"^ 29ri) = (1.4 x 10-^° eV^, 0.43, 



9.9 X 10-*' eV^, 4.0 x IQ-^) (the black triangle-up in Figs, la and 2a) Q- 

In the "up" region in the ^m^^ — sin^ 2^12 plane one has X^in — 0-03 and this minimum 
is located at (Am^^, sin2 2ei2, Am^^, sin2 20i3) = (5.8 x 10-^° eV^, 0.43, 9.9 x IQ-^ eV^, 
3.8 X 10-3). 

For E < 0.41 MeV and values of Arrig^ and sin^ 2^13 from the solution intervals (36a) 
and (36b), one has P2MSW — 1 ~ 1/2 sin^ 26*13 — 1; and consequently we obtain from eq. (9): 
P(z/e — > z/e;t£;,to) = P2voi^e "~^ i^e'-,tE-,tcd)- This implies that the effect of the MSW transi- 
tions for the pp neutrinos is negligible and they effectively take part in vacuum oscillations 
only. Moreover, for E < 0.41 MeV and Am^^ ^ 3 x 10"" eV^ we find L^^ ^ 3.4 x lO'^ km 
and 27ri?/L^i ^ 27. Correspondingly, for E < 0.41 MeV the term cos{27rR/L^^) in P^yl is a 
fastly oscillating function of R/p, which is averaged practically to zero, e.g., by the integra- 
tion over the pp neutrino energy, or by averaging over an uncertainty AE in the measured 
value of the energy E, AE/E 3> 4 x 10"^. Therefore for Am|^ ^ 3 x 10"^^ eV'^ the pp 
neutrino contribution to the event rate of the Ga-Ge detectors, for instance, is suppressed 
by the energy-independent factor P{ue — > Ue] ^e, ^0) — AVo — 1 ^ 1/2 sin^ 26*12. 

The 0.862 MeV ^Be neutrinos take part in adiabatic MSW transitions in the Sun, while 
the ^B neutrinos with E ^ 4 MeV undergo nonadiabatic transitions. Both the ^Be and ^B 
neutrinos, as well as the z/^ and/or Ur into which a fraction of the z/g has been converted by 
the MSW effect in the Sun, participate in vacuum oscillations after leaving the Sun. These 
oscillations are modulated by the MSW probability -^2^51^ (^^^ Figs. 4a - 4b). Let us note 
that for Am2i ^ 3.0 x 10~^° eV"^ the effects of the averaging of P^i^e -^ '^e]tE,to) over the 
time period of 1 year is negligible for the 0.862 MeV ^Be neutrinos; this effect becomes 
important for the ^Be z/g for Am2i ^ 5.0 x 10"^° eV"^, as Fig. 4b illustrates. For the range 
of Am2i of interest, Am'^i ^ (5.0 — 6.0) x 10"^'^ eV"^, the time averaging does not change the 



^^It is very close to the lowest x^ value possible, Xmin — 0.024, with the set of the data [||-^ we 
use in the present analysis. 
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probability P(z/e -^ ^e'-, ts, to) for solar neutrino energies E ^ 3.0 MeV, i.e., for the dominant 
fraction of the ^B neutrinos. 



With respect to the predictions in ref. |39|, the signal in the Kamiokande detector and 



the contribution of the ^B neutrinos to the signals in the Cl-Ar detector are smaller typically 
by factors of ~ (0.43 - 0.47) and ~ (0.32 - 0.36). The pp and the 0.862 MeV ^Be u^ fluxes 
are suppressed by factors of ~ (0.65 — 0.90) and ~ (0.11 — 0.27) for most of the values of 
the parameters from the allowed region. However, for sin^ 29i2 ~ 0.9, for instance, one has 
Pl^e -^ i>e;tE,to) — 0.55 for the pp neutrinos. Even in this case the 0.862 MeV ^Be Ue flux 
is reduced by a factor of ~ 0.3, but the indicated possibility is rather marginal. 

There are rather large and distinctive distortions of the ^B u^. spectrum (see Figs. 4a and 
4b) in the case of solution B. The seasonal variations due to the vacuum oscillations of the 
signals in the Super-Kamiokande, SNO and ICARUS detectors are estimated to be smaller 
than the variations in the case of the 2z/ vacuum oscillation solution 0, except possibly in 
the small region of the Atti^i — sin^ 29i2 plane where Am2i ^ 10~^° eV'^ and sin^ 29i2 ^ 0.7. 
The range of the predicted values of the day-night asymmetry in these detectors is different 
from the one expected for the 2z/ MSW solution. The seasonal variation of the 0.862 MeV 
^Be z/g flux caused by the vacuum oscillations is expected to be considerably smaller than in 
the 2z/ case, while the day-night asymmetry is estimated to be somewhat smaller than the 
one predicted for the 2z/ MSW nonadiabatic solution. The seasonal variation, nevertheless, 
may be observable. Obviously, the experimental detection both of a deviation from the 
standard (geometrical) 6.68% seasonal variation of the solar neutrino flux and of a nonzero 
day-night effect will be a proof that solar neutrinos take part in MSW transitions and vacuum 
oscillations. 

Solution C. The values of the parameters corresponding to solution C (see Figs, la and 



^^These variations were shown [^ to be not larger than 15% for the 2v vacuum oscillation 
solution. 
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2a) form a rather large region in the Arri^i — sin^ 2^13 plane, 

1.1 X 10"^ eV^ ^ Aml^ ^ 1.2 x 10"^ eV^, (38a) 

1.7 X 10"^ ^ sin2 2^i3 ^ 2.0 x 10-^ (386) 

and a relatively small one in the Arrigi — sin^ 2^12 plane, 

1.1 X 10"^° eV^ ^ Am^i ^ 1.3 x 10"^° eV^, (39a) 

0.7^ sin2 2^i2 < 1.0. (396) 

The Xmin fo^ ^his solution is larger than for solutions A and B: Xmin — 1-49 at {Am2i, 
sin2 2^i2, Aml^, sin^ 2^13) ^ (1.2 x 10"^° eV^, 0.78, 4.6 x 10'^ eV^, 5.9 x 10'^) (the black 
square in Figs, la and 2a). 

For the values (38a) and (38b) of Aml^ and sin^ 2^13 and the energies of the pp neutrinos 
of interest, we have P2MSW — 1 ^ 1/2 sin^ 2^13 = 1 and the pp neutrino flux is effectively 
suppressed by the energy-independent factor P(z/e -^ i^e',tE,to) = P2V0 — 1 — l/2sin^26'i2 
(Figs. 5a and 5b). The energy of the dominant fraction of ^Be neutrinos, 0.862 MeV, is 
for this solution in the region of energies for which the effect of the time-averaging of the 
probability P{i'e -^ '^e]tE,to) is important. This effect is negligible for E ^ 2 MeV, i.e., 
for most of the ^B neutrinos. The 0.862 MeV ^Be z/g flux is suppressed either due to the 
vacuum oscillations only (Fig. 5b) or due to the combined effect of the vacuum oscillations 
and MSW transitions (Fig. 5a). The suppression of the ^B z/g flux is a result of the interplay 
of the vacuum oscillations and the MSW transitions (Figs. 5a and 5b). As a consequence, 
i) the *B z/g spectrum is strongly deformed and ii) one expects observably large seasonal 
variations due to the vacuum oscillations of the ^B and/or, for Amli ^ 3 x 10~^ eV"^, of the 
^Be Ue fluxes. The day-night asymmetry in the signals of the Super-Kamiokande, SNO and 
ICARUS experiments is not expected to be larger than a few percent. In the case of the ^Be 
z/g flux the asymmetry is either absent (if Am|]^ ^ 3 x 10^^ eV'^) or does not exceed a few 
percent. 
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In comparison with the predictions of ref. p^, the Kamiokande signal and the signals 



produced by the ^B neutrinos in the Cl-Ar detector are reduced respectively by factors of 
~ (0.35 — 0.41) and ~ (0.24 — 0.32), while the ranges of suppression of the pp and of the 
0.862 MeV ^Be Ue fluxes read respectively 0.50 - 0.65 and 0.17 - 0.30. 

Solution D. In the Arrig^ — sin^ 2^13 plane the region of solution D (see Figs, lb and 
2a) can formally be obtained from the one of solution A by moving the latter as a whole to 
larger values of Am^-j^ and smaller values of sin^ 26*13 and by truncating two relatively small 
sub-regions located at sin^ 26^13 ^ 0.2 and at Am|i ^ 1.1 x 10~^ eV"^. In the Am^i — sin^ 2^12 
plane the solution D region is determined by 

3.2 X 10^^^ eV^ ^ Aml-^ ^ 4.0 x 10"^^ e1/^ (40a) 

0.62 ^ sin^ 2^12 < 1.0. (406) 

The minimum x^ value for this solution is similar to that for solution C, Xmin — 1-35; 
it is located at (Am^^, sin2 2^i2, Am^^, sin2 20i3) = (3.4 x lO^^^ eV^, 0.80, 1.5 x 10"^ eV^, 
2.5 X 10"^) (the black diamond in Figs, lb and 2a). 

For the values of Am^^ given in eq. (40a) and 0.22 MeV <E< 0.41 MeV, the effect of 
the time averaging of the probability Piye — ^ ^'e; ^e, to) cannot be neglected; for E ^ 0.7 MeV 
it becomes inessential. The Ue survival probability for the 0.862 MeV ^Be neutrinos is given 
by expression (35). The same expression is approximately valid for the pp neutrinos (see 
Figs. 6a - 6c). However, the pp neutrino contribution to the signals in the Ga-Ge detectors is 
suppressed by the factor given in eq. (21). The ^B z/e flux is suppressed due to the combined 
effect of the MSW transitions and vacuum oscillations (Figs. 6a - 6c). 

With respect to the predictions in ref. I^S], the contributions of the pp and ^Be (*B) 
neutrinos in the signals of the Ga-Ge (Cl-Ar) detectors are smaller by factors of ~ (0.50 — 
0.70) and ~ (0.31 - 0.49) (~ (0.22 - 0.28)), while the signal in the Kamiokande detector is 
reduced by a factor of ~ (0.32 - 040). 

As in the case of solutions B and C, the ^B Ue spectrum is strongly deformed by the 
combined effect of the MSW transitions and the vacuum oscillations. The seasonal variations 
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of the pp and of the 0.862 MeV ^Be z/g fluxes generated by the vacuum oscillations are 
estimated to be greater than ~ 10%; they are expected to be smaller than a few percent 
for the ^B z/g having energy E > 5 MeV. The day-night asymmetry in the signals of 
the Super-Kamiokande, SNO and ICARUS detectors may be observable only for values of 
sin^20i3 ^ 10^^; even in this region it is not expected to be bigger than a few percent. The 
day-night effect for the 0.862 MeV ^Be u^ flux is negligible. 

Solution D has as a ^13 ^ limit the second "low" ®B z/g flux 2z/ vacuum oscillation 
solution discussed in ref. Jill: the regions of values of the vacuum oscillation parameters of 
the two solutions practically coincide. This 2z/ solution was found to be possible for values 
of the initial ^B z/g flux which are by a factor ~ (0.45 — 0.65) (~ (0.39 — 0.56)) smaller than 
the flux predicted in ref. [^ (in ref. [^). The value of xLin fo^ ^^e indicated purely 2z/ 
vacuum oscillation solution, xLin — 5-0 (2 d.f.) [|lT|, is considerably larger than Xmin ^^ ^^e 
solution D. The MSW + VO effects and correspondingly the purely VO effects on the ^Be 
and/or ^B neutrino fluxes in the cases of the two solutions are also very different. 

Solution E. This solution holds for relatively large values of sin^ 2^12 and sin^ 2^13, 

0.70 ^ sin2 20i2 ^ 0.80, (41a) 

0.14^ sin2 2^i3 ^0.39, (416) 

and values of Ar/igi forming a tiny region around the point Am2i — 9.4 x 10~^° eV"^, with 
the values of Am|]^ belonging to the narrow interval 1.2 x 10~^ eV"^ ^ ^"^31 ^ 1-4 x 10~^ eV"^ 
(Figs, lb and 2b). The xLin value is rather large, Xmin — 3-3, and we shall not discuss this 
solution in greater detail. Let us mention only that the pp and ^Be z/g fluxes are suppressed 
by one and the same energy-independent factor given in eq. (21) - the average z^e survival 
probability in the case of 3z/ vacuum oscillations, while both MSW transitions and vacuum 
oscillations contribute to the suppression of the ^B z/g flux. This is illustrated in Fig. 7. 

Finally, we would like to point out that, as Fig. 8 indicates, the MSW transitions + 
vacuum oscillations may provide an alternative to the purely MSW mechanism of suppression 
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of the solar v^. flux in the case discussed in ref. [|^ . 
4. Conclusions. 

Assuming three flavour neutrino mixing takes place in vacuum, we have investigated the 
possibility that the observed depletion of the solar z/g flux is caused by an interplay of MSW 
transitions of the solar v^. in the Sun followed by long wave length neutrino oscillations 
in vacuum on the way from the Sun to the Earth. Choosing the neutrino mass squared 
difference A?7i|j^ > to be responsible for the MSW transitions in the Sun, the long wave 
length vacuum oscillations can be due either to Arrij]^ > (inequality (4) is valid) or to 
Am32 > ((5) holds). In both cases the solar v^. survival probability is described, as we have 
shown, by simple analytic expressions, eq. (9) and eq. (23a). However, the two cases were 
found to be very different in their physics implications. In particular, if (5) is realized, the 
vacuum oscillations are not suppressed, in principle, only for large values of the MSW mixing 
angle parameter. The case of (4) appeared to provide a much richer spectrum of possible new 
solutions of the solar neutrino problem of the hybrid MSW transitions + vacuum oscillations 
type and we have analyzed only this case. We have found for the values of Atti^i and ATTig^ 
within the ranges given in (7) and (8) several versions of this solution, denoted by A, B, C, D, 
and E in the text, most of which are distinctly different from the two- neutrino mixing purely 
vacuum oscillation or MSW solution. In particular, the values of the MSW parameters Arrig^^ 
and sin^ 26*13 of the different solutions found lie outside the regions of the 2z/ MSW solutions 
(including the regions obtained by varying the ^B and ^Be neutrino fluxes) [pr^ -[T6[|, except in 



two cases (solutions A and B) in which very small sub-regions are located within the region 
of the 1v MSW nonadiabatic solution. Two of the solutions studied, A and D, are similar 
in what regards the effective mechanism and the magnitude of suppression of the pj9 flux 



to the two "low" B Vp flux 2v vacuum oscillation solutions discussed in ||TT|. However, the 



Xmm values for the former are much smaller than for the latter. In addition, the MSW -|- 



VO solutions A and D and the two corresponding purely vacuum oscillation solutions 111 
imply very different mechanisms of suppression of the ^Be and/or the ^B neutrino fluxes. 
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A general feature of the MSW + VO solutions studied by us is that the pp z/g flux 
is suppressed (albeit not strongly - by a factor not smaller than 0.5) primarily due to the 
vacuum oscillations of the z/g, the suppression of the 0.862 MeV ^Be z/g flux is caused either 
by the vacuum oscillations or by the combined effect of the MSW transitions and the vacuum 
oscillations, while the ^B z/g flux is suppressed either due to the MSW transitions only or by 
the interplay of the MSW transitions in the Sun and the oscillations in vacuum on the way 
to the Earth. The solutions differ in the way the pp, ^Be and the *B neutrinos are affected 
by the z/g MSW transitions and/or the oscillations in vacuum. 

Searching for MSW + VO solutions we did not scan the entire region of the parameter 
space defined by eqs. (7), (8), (32) and (33). In particular, solutions which correspond 
to a suppression (by a factor ^ 0.5) of the pp neutrino flux due to the MSW transitions, 
concomitant with the requisite suppression of the ^Be and of the ^B neutrino fluxes due 
to the MSW effect and/or the vacuum oscillations and respectively the vacuum oscillations 
alone, are possible for Amj-^ ~ (10^^ - 10^^) eV^ and Am^^ ~ (0.3 - 1.0) 10"^° eV^. 

For all MSW + VO solutions we have considered, the ^B z/g spectrum is predicted to 
be rather strongly deformed. The ®B neutrinos undergo nonadiabatic transitions for values 
of the MSW parameters Am|i = (0.1 - 1.2) x 10"^ eV^ and sin^ 2^13 = (0.2 - 3.0) x lO"-'^ 
which are respectively larger and smaller than the values corresponding to the 2z/ MSW 
nonadiabatic solution. The ®B z/g adiabatic transitions take place for Arrig^ = (1.1 — 1.5) x 
10"'^ eV"^ and sin^ 2^13 = (3 x 10^^ — 0.5). Clearly, the Kamiokande data on the shape of the 
®B neutrino spectrum can be used to further constrain the solutions we have found. Such 
an analysis, however, lies outside the scope of the present study. 

For the MSW + VO solutions considered by us the day-night asymmetry in the signals 
of the detectors sensitive only to ^B or ^Be neutrinos are estimated to be rather small, not 
exceeding a few percent. The seasonal variation effect caused by the vacuum oscillations 
can be observable for ^Be neutrinos and, for certain relatively small regions of the allowed 
values of the parameters, can also be observable for the ^B or for the pp neutrinos if the pp 
neutrino flux is measured with detectors like HELLAZ or HERON. 
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Finally, the MSW transitions + vacuum oscillations can be an alternative to the purely 
MS W mechanism of suppression of the solar z/g flux in the case considered in ref . H2| . 
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Appendix. Derivation of the Expression for the Solar z/g Survival Probability 

In this Appendix we present a derivation of expression (9) for the average probability 
P(z/e —>■ Vf^^tEito), which is not based on the general analytic result for the average solar 
Ve survival probability in the case of three-neutrino mixing -P(z/e — > i^e;t,to); derived in ref. 
|2^. We remind the reader that the two relevant MSW transition and vacuum oscillation 
parameters Am|]^ and Amli (or Amig) are assumed to satisfy inequality (4) (or (5)) and to 
lie in the intervals (8) and (7), respectively. 

The probability amplitude that the solar z/e will not be converted into z/^ and/or Vj. when 
it propagates from the central part of the Sun, where it was produced at time Iq, to the 
surface of the Earth, reached at time Ie-, A^Vf. — >■ v^.^tEito), can be written in the following 
form: 

A(z/e -^ Z/g; te, to) = Yl ^Mswi^e -^ Vi; t0, tQ)Avo{vi ^ ^e\ tE, t©), (^1) 

Here t© is the time at which the neutrino reaches the surface of the Sun, AmswI^^c -^ 
z/z;t0,to) is the amplitude of the probability to find neutrino ui, I = e,fj,,r, at the surface 
of the Sun after the z/g underwent MSW transitions into z/^ and/or z/^- while propagating in 
the Sun, and Ayoiyi -^ J^e',tE,tQ) is the amplitude of the probability that ui will oscillate 
into z/g while traveling from the surface of the Sun to the surface of the Earth. For stable 
and relativistic neutrinos z/^ the amplitude Avoi^^i -^ T^e'itEito) is given by the well-known 
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expression (see, e.g., |T7| , |33| ): 



Avoii^i - ye\ tE. to) = E t/efce-^'=(*--*«)t/;, = c^^^^ ^ f/.^e-^^i^^f/;,, {A2) 

fc=i fc=i 

where i? = [Ie — t©) = Rq = 1.4966 x 10^ km, Rq being the mean Sun-Earth distance, and 
we have used Ek — Ei = Aml^/2p. 

Consider first the case specified by inequahty (4), with Am"^^ (and Am'^2) having val- 
ues in the interval (7). For p < 14.4 MeV, which is the maximal energy of the so- 
lar neutrinos detected on Earth, one has 2p/Am'^^,^2) ^ 5.6 x 10^ km, and therefore 
(Am3^(32)/2p)i?o ^ 2.7 x 10'^. As a consequence the phases {Am,l^/2p)R and (Am32/2p)i? 
in eq. (A2) give rise to fastly oscillating (with the change of p and/or R) terms in the 
probability of solar u^ survival: 

P(z/e ^ Ue;tE,to) = \A{v, ^ Ve\tEM)? ■ {A3) 

When P(z/e — > Uf,; tE, to) is averaged over (see, e.g., [0,Q) the region of neutrino production 
in the Sun, and/or the uncertainty of the energy of the neutrino detected on Earth, and/or 
over the uncertainty in the position of the solar neutrino detector [^, and/or is integrated 
over a relevant solar neutrino energy interval, the fastly oscillating terms become strongly 
suppressed and their contribution in the averaged probability, P{i'e -^ ^e), can be neglected. 
Taking this effectively into account by omitting the indicated fastly oscillating terms and 
using eqs. (Al) - (A3) one obtains an expression for the average probability, -P(z/e — > 
^e] tE, to), which is convenient to represent as a sum of two terms: 

P(z/e -> Ue] tE, to) = Si{tQ, to) + 52(ii?, t©, to) = Si + S2 , (AA) 

where 

3 

Sl = Yl J2 \Uek?{Ui^,AMSw{.^e-'l^htQ,to)){Ui,t,AMSw{.^e-> l^V]tQ,to))* , (A5) 

fc=l l,l'=e,ii,T 
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This is relevant only for the radiochemical detectors (Homestake ||l|, GALLEX p] and SAGE 
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■ 21 D 

l,l'=e,fj,,T 

{A6) 
and the bar on the quantities 6*1^2 in eq. (A4) means that all fastly oscillating terms in 
Si^2 arising from the amplitudes Amsw^^c -^ ^h^Qi to), which are rendered negligible by the 
averaging over the region of neutrino production in the Sun, etc. should be dropped. Obvi- 
ously, S'i(t0,to) is determined only by the solar z/g transitions in the Sun, while 5'2(t£;, t©, to) 
contains all the dependence of -P(z^e ^ ^e] ^Ei to) on the long wave length vacuum oscilla- 
tions between the Sun and the Earth. Note, however, that S2itE,tQ,tQ) depends on the Ue 
transitions in the Sun as well. 

The MSW transition amplitudes AMswi^e -^ ^ht,to) = Ai{t), t <tQ, I = e,fi,T, satisfy 
the system of evolution equations: 

^|^'W= E Mait)A,{t). {A7) 

l'=e,fi,T 



Here M(t) is the evolution matrix which can be chosen in the form (see, e.g., |19,20[|): 



Mii'it) = ;^E Ui,A^%U+ + A{t)6,i6a' - SwiY^ \UekfAml, + A{t))], {A8) 

-^P j=2 k=2 

where 

A{t) = 2pV2GFNe{t), {A9) 

and Ne{t) is the value of the electron number density at the point of the neutrino trajectory 
in the Sun, reached at time t. The initial conditions for the system of equations (A8) read: 

Ae(to) = 1, A^ito) = Ar{to) = 0. (^10) 

Since Arri^i < 5 x 10^^° eV"^ and Arrig^ > 10"''' eV"^, we can expect that the evolution 
of the neutrino system in the Sun practically does not depend on Am^i. At the same 
time the lower energy pp neutrinos, for instance, can take part in vacuum oscillations in 
the Sun induced by Am^i p2|,p5|,p6|. Inspecting eq. (A8) for the evolution matrix M{t) 



one could conclude that the terms containing Ar/igi ^^ ^ factor can be neglected provided 
\Ui3Aml^U;,^\ > \Ui2Aml^U;,2\, I ^ V = e,fi,T, and \{\Ui3? - \Ue3\^)Aml^\ > \{\Ui2\^ - 
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|f/e2p)Am|]^|, / = e,/i. Actually, using some of the intermediate results of the analysis 
performed in ref. [^ one can convince oneself that for Aml^ and Am|j^ having values in the 
intervals (7) and (8) the effects of Arn^i on the solar neutrino transitions/oscillations in the 
Sun are either effectively accounted for by the term S'2(tB, t©, to) or are negligible (even for 
the pp neutrinos). This implies that one can formally set Ar?7,2i to zero in the expression for 
M{t), eq. (A8). In this case the evolution matrix M{t) depends only on the elements Ui^, 
/ = e, /i, r, forming the third column of the matrix U. 

The phases of the elements U13, U13 = e'^°''-^\Ui3\, where a/3 are real constants, can be 
eliminated from M{t) by a redefinition of the amplitudes Ai{t): Ai{t) -^ e~*°'^y4/(t), / = 
e, n, T. The evolution matrix M{t) of the system of evolution equations for the amplitudes 
e~^°'^^Ai{t) can formally be obtained from M{t) by replacing Ui^ by {Ui^l in the expression 
(AS) for M{t). As we shall show, the average solar z/g survival probability -P(t'e ~^ '^e',tE,to) 
does not depend in the case of interest on the phases of the elements of the lepton mixing 
matrix U, and, in particular, on ai^. 

Performing the transformation 



V 






e-'''-^A^{t) 



( 



1 












KM. 



V|t/,3P + |C/p3P V|C/r3P + |C/M3l' 



I^/msI 



l^ral 






V 



-«ae3 



A;,(t) 



(All) 



V \J\VrZ? + \U^aV V|C/r3p + |C/M3p / 

we find that the amplitudes (Ae(t),A|^,(t),A'^,(t)) satisfy a system of evolution equations 
having the same form as the system of equations (A7) with an evolution matrix M'{t) given 
by the expression: 



^'^'W - Tp 








U,.\Jl-\U..\^ Am. 



e3 



'31 







\Ue3\^Aml^-A{t) 







f/.,U/l-|[/.,P Am^ 



e3| 



''31 







;i-2|f/e3p)Amii-A(t) 



(A12) 



In obtaining eq. (A12) we have used the unitarity of the lepton mixing matrix U: |f/e3p + 



|f^u3p + l^-r3p = 1- The initial conditions for the new system of equations follow from eqs. 
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(AlO) and (All): 

Ae(to) = l,A'^'{to) = K>{to) = 0. (A13) 

It follows from eq. (A12) that the evolution of the amplitude A'^f(t) is decoupled from 
the evolution of Af.{t) and A'^,{t), while Ae{t) and A'^i{t) satisfy a two-neutrino mixing 
system of evolution equations describing two-neutrino MSW transitions of the solar Ug in the 
Sun. The transitions can be enhanced by the solar matter effects, as it follows from (A12), 
provided |f/e3p < 0.5. Explicit expressions for the amplitudes Ae{t) and A'^i{t) have been 
obtained, in particular, assuming Ne{t) decreases exponentially along the neutrino path in 
the Sun [^,^,Q. In this case the system of evolution equations for A^it) and A'^i(t) can 
be solved exactly [^,0. Let us note that the exponential approximation describes with a 



good precision the change of Ne{t) along the neutrino path, predicted by the solar models 

in- 





A',>it) = 0, 


A,it) = 


u,. 


^f^MsP+lf^rsP 


Ar{t) = 


f/.3 


J\Ur,\^+\U,,^ 



As it is easy to see, eqs. (All) - (A13) (see also (A8)) imply: 

(A14a) 
-"'-'A'^,{t), {AUh) 

-^"=M;,(t). (A14c) 

Evidently, one has: 

\AM^ + \AM^ + |A.(t)P = \AM? + \K'{t)? = 1- (^15) 

We shall express next the probability P{ve -~* ^e',tE,tQ) in terms of the two amplitudes 
Ae{t) and A'^'{t) which contain the whole information about the MSW transitions of neutrinos 
in the Sun. Inserting expressions (A14b) and (A14c) for A^{t) and A,-(t) in equations (A5) 
and (A6) and using the unitarity of the lepton mixing matrix U we obtain very simple 
expressions for the quantities Si and S2: 

S,{tQ,to) = [(1 - \U,3\y - 2\UMUe2\'] |^(t0) - , '^'^'' A'Ate)\' 



lUesl^ 
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+ If/esr \Ae{tQ)\Ue3\ + Jl - jf/esP ^'(^0)!', (AI6) 



l^esl ./ /, x|2 ~ -R 



S2{tE,tQ,t^)=2\UMUe2? I^fe)- , ' ^'' K,{tQ)Y cos27r^^. (A17) 



Let us denote 

|Ae(t)P = 1 - |A;,(t)|2 = P2S^^(l/e ^ ^/e;t0,to) . (A18) 

Explicit expression for the probability P2Mswi^e ~^ '^eSioj^o) can be found in ref. [^. 
However, as we have already emphasized, the oscillating terms in P-^uswi^e -^ ^'eS^o^^o) 



give a negligible contribution in the averaged solar z/g survival probability ||3^ . Therefore 
only the average term P-zuswi^e — ^ i^e] te, to) in P2Mswi^e ~^ ^e'l '^o? ^o) should be taken into 
account when one obtains 5*1(^0,^0) and S2(tE,tQ,to) from eqs. (A16) - (A18). 

The last quantity to be computed is the interference term Re[Al{tQ)A'^,(tQ)] in eqs. (A16) 
and (A17). This term coincides with the quantity Re[i?H(to, to)] in eq. (29) in ref. [^] when 
tp = t©. An explicit analytic expression for Re(y4*(t0)y4!^,(t0)) can be found exploiting, e.g., 
the exact analytic results for Ae{tQ) and A'^,(tQ) derived in ref. [^ (see eqs. (28) and (29)) 
and cast in a form which is more convenient to use in calculations like the present one in 
ref. [p5l , p8| . It follows from eqs. (11), (12) and (18) - (20) in ref. ||3^, for instance, that 



Re[A*(t0)v4;,(t0)] = - (- - P('3^)) sin20i3 cos20^(to) + Oscillating terms, {A19) 

where the angles ^13, ^^(to) and the "jump" probability P^^^n are defined in eqs. (12), (13) 
and (15). The oscillating terms appearing in eq. (A19) do not contribute in the averaged 
solar Ue survival probability]^. Let us note that, as can be shown [Q, the result eq. (A19) for 
Re[A*(t0)A!^, (t0)] is general: the fact that we have used the specific (but exact) exponential 
density solutions for v4e(t0) and A'^,(tQ) to derive it is refiected only in the specific form of 
the "jump" probability P/g^-, (eq. (15)) one gets in eq. (A19). 
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Explicit analytic expression for the oscillating terms in eq. (A19) can be obtained from the 



formulae given in refs. [p^,35.38| 
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Taking into account the remarks we have made concerning the oscillating terms in eqs. 
(A18) and (A19) and utilizing eqs. (A16) - (A19) and (10) - (14) it is not difficult to get eq. 
(9) from eq. (A4). 

If the Zee model |2^ type relation (5) between Am^i, Amg^ and A'm'^2 is valid, we have 
to repeat the above analysis interchanging the indices 1 and 3 in the quantities Uij and 
in Ami- entering into eqs. (A6), (A8), (All), (A12), (A14b), (Al4c), (A16) - (A18). In 
particular, the evolution matrix M'{t) in the system of evolution equations for the amplitudes 
{Ae{t),A'^,{t),A'^>{t)) now has the form: 



It/eily/l-lt/elP Am?3 

^'W = ^ -|f/,i|2Am?3-A(t) 

If/eil^l-lf/eiP Am^g {l-2\U,,\^)Amj^-A{t) 

{A20) 

where Amf^ < 0. The new form of M'{t) and eqs. (A16) and (A17) in which |f/e3(i) P and L21 

are now replaced respectively by |f/ei(3)P and L23, entail the introduction of the new angles 

^23, ^13 and O^ito) defined by eqs. (24) - (26), as well as of the "jump" probability Pn^y 

It also leads to the requirement |f/eiP > 0.5 which is a necessary condition for the MSW 

resonance. Comparing the evolution matrices (A12) and (A20) it is easy to deduce the form 

of PL^y given the expression for -P(3i), eq. (15). The interference term Re[Al(tQ)A'^,(tQ)] is 

given by eq. (A19) in which 6*13, 0^(to) and PL^^ are replaced by 6[^, O^ito) and -Pm)- Using 

the above remarks it is easy to get expression (24a) for the average probability P^{ve -^ 

Ve\ tE,to)- 
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Figure Captions 

Figs, la — lb. Regions of values of the MSW transition parameters ^rn"^^ and sin^ 26'i3 for 
wliicli tlie x^— analysis of the solar neutrino data gives x^ ^ 3.841. The regions shown are 
obtained for fixed values of the vacuum oscillation parameters A?7i|^ and sin^ 29i2 which are 
indicated in the figures. Shown are also the 95% CL. 2u MSW nonadiabatic and adiabatic 
solution regions. 

Figs. 2a — 2b. Regions of values of the vacuum oscillation parameters A?7i2i and sin^ 29i2 
allowed by the solar neutrino data (x^ < 3.841). The regions correspond to fixed values of 
the MSW transition parameters Amg^ and sin^ 2^13 indicated in the figures. 

Figs. 3a — 3c. The averaged probability -P(z/e -^ z/g) (solid line) as a function of the 
solar neutrino energy E for four sets of values of the parameters indicated in the figures 
and belonging to region A. The dependence on E of the MSW probability P2Mswi^e ~^ ^e) 
averaged over the region of the ^B neutrino production is also shown (dashed line). 

Figs. 4a — 4b. The probability -P(z^e —>■ ^e) (solid line) not averaged (a) and averaged (b) 
over a period of 1 year |]1T| (see the text), as a function of the neutrino energy E for two 
different sets of values of the four parameters from region B. The averaged MSW probability 
P2Mswi^e — > Ue) as a fuuctiou of E is also shown (dashed line). 

Figs. 5a — 5b, 6a - 6c, 7. Examples of the dependence of the averaged probabilities 

— — ro-\ \ 

P{h'e —^ i^e) (solid line) and P-iuswi^e — ^ t'e) (dashed line) on the neutrino energy E in the 
cases of solutions C, D and E, respectively, considered in the text. 

Fig. 8. The averaged probability -P(z^e —^ '^e) (solid line) as a function of E for values of the 
parameters (indicated in the figure), for which the MSW transitions + vacuum oscillations 
can provide an alternative to the purely MSW mechanism of suppression of the solar i/g flux 
in the case discussed in ref. 1021. 
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